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Abstract 

We discuss various ways in which the computation of conservative 
Gravitational Self Force (GSF) effects on a point mass moving in a Schwarz- 
schild background can inform us about the basic building blocks of the Ef- 
fective One-Body (BOB) Ifamiltonian. We display the information which 
can be extracted from the recently published GSF calculation of the first- 
GSF-order shift of the orbital frequency of the last stable circular orbit, 
and we combine this information with the one recently obtained by com- 
paring the EOB formalism to high-accuracy numerical relativity (NR) 
data on coalescing binary black holes. The information coming from GSF 
data helps to break the degeneracy (among some EOB parameters) which 
was left after using comparable-mass NR data to constrain the EOB for- 
malism. We suggest various ways of obtaining more information from GSF 
computations: either by studying eccentric orbits, or by focussing on a 
special zero-binding zoom- whirl orbit. We show that logarithmic terms 
start entering the post-Newtonian expansions of various (EOB and GSF) 
functions at the fourth post-Newtonian (4PN) level, and we analytically 
compute the first logarithm entering a certain, gauge-invariant "redshift" 
GSF function (defined along the sequence of circular orbits). 



1 Introduction 



The detection of gravitational waves from coalescing binary systems depends 
upon the prior knowledge of accurate theoretical models of the emitted gravita- 
tional waveforms, so as to be able to extract the gravitational wave signal from 
the noisy output of the detector. There has been much progress, over the past 
few years, on the development of accurate computational tools for describing the 
motion and radiation of (comparable-mass) compact binary systems (i.e. sys- 
tems made of black holes or neutron stars). These computational tools are based 
either on analytical methods, on numerical ones, or on various combinations of 
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both. The first formahsm which made several quantitative and quahtative pre- 
dictions about the entire coalescence process of comparable-mass circularized 
black hole binaries (from early inspiral to late ringing) is the analytical Effec- 
tive One Body (EOB) formalism [II HI [3]. [Note that the EOB formalism uses, 
as essential inputs, the results of high-order post-Newtonian (PN) expanded 
results (see [1] for a review). However, it does not use PN results in their 
original "Taylor-expanded form" , but rather in some suitably resummed form.] 
Soon afterwards, a combination of (short) full numerical simulations, with a 
"close limit approximation" [5] to the ringing final black hole led to the first, 
numerical-based, (approximate) description of the coalescence of (comparable 
mass) circularized black hole binaries [51 [7] . Recently, several breakthroughs in 
numerical relativity (NR) [8] have allowed numerical methods to describe, with 
very high accuracy, the motion and radiation of coalescing black holes. These 
impressive NR achievements do not, however, render obsolete the development 
of analytical methods for describing the motion and radiation of coalescing black 
holes. Indeed, in spite of the high computer power used in NR simulations, the 
calculation of one sufficiently long waveform (corresponding to specific values 
of the many continuous parameters describing the considered binary system) 
takes on the order of several weeks. For detection purposes, one needs to com- 
pute tens of thousands of theoretical templates, so as to densely sample the full 
parameter space. This is a clear motivation for developing accurate analytical 
models of waveforms. One avenue for doing so is to use the natural flexibil- 
ity of the EOB formalism (which was emphasized early on [3]) to "tune" some 
of the theoretical EOB parameters representing yet uncalculated, higher-order 
effects until the EOB waveform "best fits" a sparse sample of high-accuracy 
NR waveforms. Over the last years, this strategy has been vigorously pursued 
and has led to an impressive analytical/numerical agreement, with residual dif- 
ferences in phase and amplitude on the order of the current numerical errors 
[9l [TOl [TTl [37l [131 [Ml [15] . The most recent, and most accurate, implementation 
of this strategy used only two EOB flexibility parameters, denoted 05 and Og, 
and found that there is a strong degeneracy between 05 and ag in the sense that 
one can find an excellent NR-EOB agreement within a long and thin banana- 
like region in the (05,05) plane [T31[TB]. This "good fit region" approximately 
extends between the points (05,05) — (0,-20) and (05,05) = (—36, +520) [16] . 

The main purpose of the present work is to study to what extent the com- 
putation of gravitational self force (GSF) effects on a point mass moving in a 
Schwarzschild background can inform us about some of the yet uncalculated 
higher-order theoretical EOB parameters, such as the just mentionned 05 and 
og parameters. The GSF program aims at describing the motion and radia- 
tion of a small compact object (of mass mi) moving around (and eventually 
inspiralling into) a much larger central black hole (of mass 7712). By contrast to 
the computational tools mentionned above (PN, EOB, NR), which can study 
binary systems with arbitrary (symmetric) mass raticll] v = TO1TO2/ {"mi + "t-2)^, 

^Note, however, that NR simulations become increasingly difficult and time consuming as 
the symmetric mass ratio v = mim2/{m\ + m2)'^ becomes small. 
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the GSF program is a priori limited to the extreme mass ratio case, mi <^ TO2, 
i.e. 1/ <C 1. Roughly speaking, the GSF program is motivated by the planned 
space- based interferometric gravitational wave detector LISA [T7], while the 
other computational programs are motivated by the higher-frequency ground- 
based interferometric detectors (such as LIGO |IH])- The GSF program has 
been under development for a long time (essentially since the classic work of 
Refs. [m [20l [m [22l [23]), but has started to reach fruition only very recently. 
[See for recent reviews of the GSF program.] We have particularly in 

mind here the recent breakthrough of Barack and Sago (BS), Ref. which 
computed the first-order GSF correction to the frequency of the Last (circu- 
lar) Stable Orbit (LSO) of a small point mass mi orbiting around a large (non 
spinning) black hole of mass m2. [Another interesting recent result is the com- 
putation of the first-order GSF correction to a certain gauge invariant function 
of the sequence of circular orbits [27l [28j .] We wish here to relate the published 
result of [26] to the parameters entering the FOB formalism (and notably the 
parameters 05, ag, . . . mentionned above), and suggest further ways of using the 
GSF program for learning more about the theoretical building blocks of the FOB 
formalism. Hopefully, such a bringing together of the GSF and FOB programs 
might benefit both, and might thereby help one to develop better theoretical 
models of coalescing binaries. Note finally that we shall only consider here the 
conservative aspects of both the EOB and the GSF. Indeed, the result of [55] 
only concerns the time-symmetric, conservative piece of the GSF. Note that the 
influence, due to conservative interactions, of a non vanishing symmetric mass 
ratio ly (with < < 1 /4) on the frequency of the LSO was first studied within 
the EOB formalism in [1] at the 2PN level, and in [29] at the 3PN level. As 
for the influence of time-odd, radiative interactions on the location of the LSO 
it was first studied in [5] (for arbitrary 1/ < 1/4) and shown there to lead to a 
"blurring" of the LSO frequency of fractional order v'^^^. [See also the work of 
[30] . restricted to the small mass ratio case <C 1.] 



2 Notation, choice of coordinate system, and 
physical units 

Let us start by warning the reader about a conflict of notation, between the 
standard EOB notation and a usual GSF one, which can bring confusion. We 
shall here adhere to the standard EOB notation in which the two masses of the 
binary system are denoted mi and 1112 (with, say, mi < m2 to fix ideas) and 
where one then defines 



M = mi+m2; fJ. = ■ ; ly = — ^ ■ . 2.1) 

1711+1712 M (mi+m2) 



By contrast, many GSF works use the letter M to denote the large mass (de- 
noted m2 in our notation), and the letter to denote the small mass (mi in our 
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notationjl. 

Another possible source of confusion concerns the choice of coordinate sys- 
tem to describe the binary system. Though we shall only work with (formally) 
"gauge invariant" quantities, this still leaves some potentially confusing ambi- 
guity. Indeed, in PN and EOB (as well as NR) studies, it is tacitly assumed that 
one uses an "asymptotically flat coordinate system" , i.e. a coordinate system 
which exhibits, in a standard mannei[f|, the asymptotic flatness (g^^ — > rjfj^^) 
of the metric generated by the considered binary system. By contrast, while 
some GSF works (e.g. that of [23 based on a "Regge- Wheeler gauge") do use 
"asymptotically flat" coordinates, the GSF works that use a "Lorenz gauge" 
for studying perturbations h^^ix^) off a Schwarzschild background g^^}{x'^) (as 
is the case of the work ^26^ that we shall consider) actually use a coordinate 
system which is not (explicitly) "asymptotically flat". Indeed, it was shown 
by Detweiler and Poisson 31], and by Barack and Lousto |32| (building up on 
results of Zerilli [22,) that the unique low multipole {(. < 2) contributions to 
the first-order Lorenz gauge metric perturbation generated by a point mass in 
circular motion are such that the monopole {£ = 0) contribution (and only the 
monopole) yields a metric perturbation which does not decay as r ^ oo, but 
tends to 

\imh]:r'^^ix^)dx^dx'^ = -2 (2.2) 
r-^oo ro(l - 2Gm2/ro) 

+ dr^ +r^de'^ + r^sm^ed(p'^ + 0{mf). 

Here, rg is the Schwarzschild (areal) radial coordinate of the considered circular 
orbit, and 

E,^^^-ul, (2.3) 
mi 

the (toi) particle's conserved energy per unit rest mass. Note that we generally 
use units such that c—1 (and often also G = 1), and a "mostly plus" signature 
( h ++). Adding the Lorenz-gauge perturbation (|2.2p to the ("asymptoti- 
cally flat" ) background Schwarzschild metric leads to a perturbed metric for 
the binary system such that 

lim(g(°jSchw^/^Lorcnz)^^^^^. ^ -(l + 2a)diL,_ (2.4) 

+ dr"^ +r'^de'^ ^-r"^ sir? edLp^ + 0{ml) , 



where 



Gmi E, ^2.5) 



ro(l - 2Gm2/ro) 



^We find helpful to use the mnemonics 1 < 2 to remember that, by dofiniton, mi < r?i2 
(and mi <C m2 in the extreme mass ratio case). 

^Possibly after a trivial transformation from Minkowskian-like to other coordinates, say of 
the (t, r, 8, ip) type. 
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In the case that we consider here of circular orbits (and that will suffice for our 
purpose), the explicit expression, in terms of the orbital radius tq, of the specific 
conserved energy Ei is 

E, = ] , (2.6) 

so that (|2.5p can be reexpressed as 

Gm^ 

a= (2.7) 

ro^yi- 3Gm2/ro 

The result (|2.4p is not new, and played, in particular, a crucial role in a recent 
work |33j which compared the numerical results of the GSF obtained by two 
different methods, using different gauges. We are spelling it out again explicitly 
here to emphasize the following point. 

Through first order in mi , we need to renormalize the time coordinate asso- 
ciated to Lorenz-gauge GSF calculations to work in normal, "asymptotically flat 
coordinates" . Namely, in the case considered here of a point mass mi moving 
on a circular orbit of radius rp, we need to introduce 

ifiat = (1 + a + C'(m^))tL orcnz , V 

together with a^^at = r sin 6* cos (/?, ?/flat — r-sin6'siniy9, 2;flat = rcosd, to ensure 
that the perturbed metric ds^ — {gjS} + h^j^^) dx'^dx" tends, as r ^ cxd, towards 
the standard Poincare-Minkowski metric lim ds^ = —dti ^ + dxi . + dy'i ^ + 

r — >oo 

rf^flaf ^ standard practice in relativistic gravity to use such "asymptot- 

ically fiat coordinates" because they relate coordinates to physical units. For 
instance, the asymptotically measured orbital frequency of, say, a circular or- 
bit, expressed in the physical time units (say the SI second) defined by the 
metric ds^ , Ophys, coincides with the coordinate angular frequency defined in 
"flat" coordinates fiflat = dip/dtnax- For the purpose of gravitational wave ob- 
servations, we are interested in the value of such physical-units frequencies: 
^phys = ^^flat = dip/dtfii^t. In conclusion, the physical-unit value of a particular 
frequency fiphys — f^flat differs from the corresponding Lorenz-gauge frequency 
f^Lorenz = dip / dti^oicnz by the "rcnormalization" factor deduced from 



f^phys = = 2^^ = (!-« + 0{ml)) f^Lorcnz • (2.9) 



Applying Eq. (|2.9p to the recent, Lorenz-gauge GSF computation of the LSO 
frequency [26], namely 



Gm2 S ^Lorcnz — 

with 



Rs ni-i „ / / mi 
1 + cH^ — + O — 
m2 \ \m2 , 



(2.10) 

cH^ = 0.4870(6) , (2.11) 
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(where the number in parenthesis indicates the error on the last digits), we 
conclude that the physical-unit LSO frequency is 

where the value of a, Eq. (|2.7p . at the (zeroth order) LSO, i.e. tq = 

6Gto2 + 0{mi), so that 



a 



LSO 



\/2 mi 1 mi 
6 m2 ^/T8 "12 



In other words 
Gm2 1^1;,^° =6-3/2 



1 \ TOl 



18/ "12 




(2.13) 



(2.14) 



Note that the numerical correction l/vTS ~ 0.235702 leads to a quite 
significant decrease (by nearly 50%) of the raw Lorenz-gauge coefficient c^^, 
Eq. pUll : 

1 



„BS 



18 



0.2513(6). 



(2.15) 



We recommend that, in the future, the results of any Lorenz-gauge (dimen- 
sionful) quantity be explicitly reexpressed in physical units, so as to avoid the 
need of doing such a posteriori renormalizations. A simple way of doing so 
would be to explicitly change the time coordinate, a la Eq. (|2.8p . to work in a 
"renormalized Lorenz gauge" where the metric becomes asymptotically flat. 

To compare the renormalized GSF result (|2.14p to the theoretical predic- 
tions of the EOB formalism, it is useful to introduce another "renormalization" . 
Indeed, the natural quantity for adimensionalizing a frequency in the EOB for- 
malism is {GM)~^ = {G{mi + 7712))-^, rather than (Gm2)~^, as used in GSF 
studies. This leads us to reexpressing (|2.14p as 



GA/rfLSO =G(mi 



-t- m2) f^LSO ^ 6-3/2 [1 ^ ^rcn ^ ^ ^(^2)] ^ 



(2.16) 



where the "EOB-renormalized" frequency sensitivity coefficient cfj" is given by 

1 



= 1 



= 1.2513(6) . 



(2.17) 



As a final reexpression, it is useful to replace any adimensionalized orbital 
frequency MD, by the quantity 



fGMn'^ 2/3 



(2.18) 



which plays a crucial role in PN and EOB developments. [Here, we have inserted 
the velocity of light c (elsewhere set to one) as a convenient mnemonics to 
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remember the exponent 2/3. Indeed, x ~ /c^ is the basic PN expansion 
parameter, and the exponent 2/3 converts the factor into a 1/c^ one.] In 
terms of a;, Eq. ((^1^ . Eqs. friB)) . (PT7|) become 

^LSO ^ 1 ^ ^rcn ^ ^ 0(^2)] ^ (2.19) 


with 




= 0.8342(4) . (2.20) 



It is the number c^''", Eq. (|2.20p . that we shall primarily use in the following 
to compare the GSF results to the predictions of the EOB formalism. Let 
us note in passing that all the EOB estimates [H [29] of the mass-ratio (i.e. 
v) dependence of the total- mass-adimensionalized LSO frequency GAI Vl^^^ = 
j^LSO (j^-j ^ .^^gjj j-|-g Q^j^gj. possible PN estimates [551 IM] j agree in predicting 
that 0'^^*-'(zy) is an increasing function of u, i.e. that the sensitivity coefficients 
Cq^, in Eq. (|2.16p . or (f°^, in Eq. (|2.19p . are strictly positive. We just wish to 
remark here that this property is not necessarily linked to the Lorenz-gauge, 
m2-adimensionalized sensitivity parameter c^^, Eq. (j2.10p . being positive. The 
positivity of Cj™ is equivalent to the much weaker condition c|^^ > -y= — 1 ~ 
-0.76430. 



3 A short review of the conservative EOB for- 
maUsm 

The EOB formalism comprises three main building blocks: (i) a conservative 
dynamics; (ii) the inclusion of radiation reaction effects; and (iii) the construc- 
tion of a resummed waveform. Let us review here the item (i), which suffices 
for the comparison to the corresponding conservative GSF effects. 

The conservative EOB dynamics is defined by a Hamiltonian of the form 

H^OB = M^l + 2u{H,ii - 1) , (3.1) 

where -ffcff = ^off / M denotes the "effective Hamiltonian" (per unit /i-mass) 
of an "effective particle" of mass /i, following (modulo quartic, and possibly 
higher-order, terms in the radial momentum p^), a geodesic in the "effective 
metric"0 

dslf^ = gf^{x)dx''dx'' = -A{r:v) dt^ + B{r,iy) dr^ +r^{d9^ + 81-0^ 9dip^) . (3.2) 

More precisely, the conserved energy of the effective dynamics for fi, Hcs = —po, 
is obtained by solving an "effective one-body Hamilton- Jacobi" equation of the 
form 

= fi^ + g^si^)Pf^P^ + Qip)^ (3-3) 

*The radial metric coefficient is here denoted B{r), instead of the notation B{r) used 
in the original EOB articles, because we wish to keep the notation B for another use; see 
below. 
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where — dS{x) / dx^ {S denoting the action), and where 



Q{p)^A>^''p^x)p,p,PpP^ + ... (3.4) 

denotes quartic-in- momenta, and possibly higher-order in p, contributions. As 
shown in one can, at the 3PN level, restrict the quartic correction p.4p to 
depend only on the spatial components of the momentum, and, more precisely, 
to be proportional to the fourth power of the radial momentum: Q oc p'^. We 
shall assume here that the fact that Q depends only on pr and r, and vanishes 
(at least) like pf when pr remains true at higher PN levels. Solving the 
Hamilton-Jacobi equation ()3.3p with respect to Hcs ~ —po > leads to 



iJcff = ^Air; (a*^ + J + + Q{r,Pr)) (3.5) 

where J denotes the conserved effective angular momentum. For motions in 
the equatorial plane 9 — ^ oi the coordinate system of (|3.2p . one has simply 
J — Pip. One of the basic principles of the EOB formalism is that the effective 
angular momentum J is identified with the total angular momentum of the 
binary system. By contrast, the total energy of the binary system, £ = M+ 
binding energy, is not equal to the effective energy, say, £cS — HcS, Eq. (j3.5|) . 
but to the EOB Hamiltonian (|3.ip : £ — Heob- Still, we see from p.5p that 
-ffoff governs, for a given value of J, the radial motion: 



72 2 

i/2g(r,p.; J) = A{r- v) [ + — + -^p— + Q(r,p,,) ) , (3.6) 



where, at the 3PN level. 



Q{r,pr) = 2{A-?,v)v&^^. (3.7) 

In practical calculations, it is quite useful to rescale the various quantities, 
r,pr, J,£cS into dimensionless ones, using 

t~GMi, r = GMr, pr = fipr , J=GAIfij, £cs = fi£cS- (3.8) 

Note that the effective action 

5= -£efft + JV + ^rad(r), (3.9) 

is rescaled into 5 = 5/ GM fi, so that, for instance, Pr is canonically conjugate to 
f with respect to the (rescaled) symplectic structure defined by 5. The rescaled 
version of Eqs. dM]), (IS21) then reads (at the 3PN level) 

H!s{f,Pr;j)=A{f;,.)(l + -^ + ^^ + Qif,pr)) , (3.10) 
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where, at the 3PN level, 

Q{f,Pr)^2{A-iv)v^. (3.11) 

In the present work we shall, like Ref. [26 , mostly focus on small-eccentricity 
orbits, i.e. on the limits where <C j"^ jv^ ■ Then, working to first-order in the 
squared eccentricity ~ we can neglect the quartic (or more) contribution Q 
which is ©(e^). We then see that, in this case, the conservative EOB dynamics 
depends on the knowledge of two functions of f and namely the metric func- 
tions A(r;v) and B(r;v) entering the effective metric (13.21) . In previous EOB 
works, it was found convenient to replace the second metric function B = 
by the combination 

D(r) = A{r) B{r) . (3.12) 

Here, we shall find it even more convenient to use instead the inverse of the D 
function (which naturally entered Ref. 3J, where it was proposed to resum D 
by working with the Taylor expansion of its inverse). We shall denote it as 

D{r) = D-^{r) = {A{r)B{r))-K (3.13) 

Another useful EOB notation is to introduce a special notation for the inverse 
of the scaled EOB radius f, namely 

1 GM 

u=- = . (3.14) 

r r 

With this notation in hand, the conservative EOB dynamics of small-eccentricity 
orbits depends on the knowledge of two functions: A{u\ v) and D{u\ v). At the 
3PN level, these two functions have been found to take the values [29] 

A{u]v) = l-2u + 2vu^ + vaiu'^ + 0{u^) , (3.15) 

D{u;v) = 1 + + 2{2Q-?,v)vu^ + 0{u^) , (3.16) 

where 04 = 94/3 - 4l7rV32. 

As summarized in Eqs. (|3.15p . (|3.16[) the PN calculations give access to the 
first few terms of the expansions of A and D in powers of u, for a fixed value 
of the symmetric mass ratio v. If we formally extend the definition of the EOB 
formalism, and of the A and D functions, to an arbitrary PN order (as was 
shown to be possible in Section III of [IS]), assuming that the "correction" term 
Q is kept of order as pr — > 0, we can think of A{u; v) and D{u] v) as being 
some yet unknown mathematical functions of two variables. 

At present, one has in hands two different computational tools for acquir- 
ing some knowledge of the unknown "abstract" EOB potentials A{u\ v) and 
D{u]v). On the one hand, existing PN calculations (at the 3PN level) have 
given us access to the u-expansions (|3.15p . (|3.16p (and to the first term (|3.1ip 
in the Q term). On the other hand, NR simulations gives us access (when v is 
not too different from 1/4) to various data that depend on A and D and can 
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therefore be used, in principle, to map the two-dimensional profiles of A{u; v) 
and D{u] v). The first attempt to do so dates from 2002. Ref. [35] coinpared NR 
data on the gravitational binding energy of (corotating and wavelesqj) circular 
binary systems to the EOB predictions, and looked for the best fit to a gen- 
eralized A function (|3.15p . extended by an additional (4PN) term -\-a^{v)u^. 
[Note in passing that their result (for v — 1/4) suggested a negative value 
for the coefficient of w^, a^ij) — —3, i.e. a slope as ~ —12 if one consid- 
ers a linear dependence on : 05 {v) = 05 v. Such a negative value is not 
unreasonable, as we shall further discuss below.] More recently, many works 
[ini [ni [371 m [H [IS] have compared NR data on the waveform emitted by 
radiation-reaction-driven, inspiralling binary black holes to the predictions of 
several extended versions of the EOB formalism. These extensions consisted in 
adding parametrized extra contributions not only to the PN expansion of the A 
function ((^A^'^ = v a^u^ + v a^u^), but also to the other building blocks of the 
EOB formalism: the radiation reaction force, and the emitted waveform. We 
shall compare below the results of the most recent NR/EOB comparison of this 
type, and the recent GSF result ^X^., ^T^. 

The main purpose of the present work is to explicate how GSF results can 
also be used as a third computational tool for acquiring some knowledge about 
the conservative EOB potentials A{u\v) and D{u\v). The first point to em- 
phasize is that current GSF technology can, at most, give us access to the first 
terms in the expansions of the functions Aiu;^) and D{u;v) in powers of the 
symmetric mass ratio v. More precisely, let us consider the expansions of the 
(unknown) exact functions A{u; and D{u; i') in powers of v, say 

A{u;i^) = l-2u + iya{u) + iy^a2{u) + 0{iy^), (3.17) 

5(u;i/) = 1 + 1^ d{u) + ly'^ d2{u) + 0(1^^) , (3.18) 

where, for notational simplicity, we have suppressed the index 1 on the contribu- 
tions which are linear in v. vai^u) = i'a{u), i'di{u) = vd{u). In the following, 
we shall show to what extent GSF studies can inform us about the two v-lineav 
functions of u, a{u) and d{u). 

Before plunging into the details of how GSF studies can tell us something 
about the two functions a{u) and d{u) entering the conservative EOB dynamic^j, 
let us emphasize the complementarity of the expansions (|3.15p , (|3.16p versus the 
expansions (|3.17p . (|3.18p . The PN expansions (|3.15p . (|3.16p proceed in powers 
of M = GM /c^r, but compute, at a given order in u, the exact dependence on v. 
On the other hand, the "GSF" expansions (|3.17p , (|3.18p proceed in powers of v, 
and can in principle give access, at a given order in v, to the exact dependence 
on u. For instance, the first-order GSF coefficients a{u) and d(u), if known 
exactly, would give us information about some arbitrarily high PN contributions 

^i.e. admitting a helical Killing vector. 

^To which should be added another function (of u and pr) parametrizing the linear-in-f 
piece in the the "exact" version of the Q contribution in Eq. 113.101 1. As said above, we shall 
consider here observables that do not depend on the Q contribution. 
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in (|3.15p . (|3.16p . Reciprocally, as the PN results p.l5p . (|3.16p have made no 
truncation on the powers of i^, we see that Eqs. p.lSp . (|3.16p is already giving us 
information about high-order terms in the "GSF expansions" (|3.17p . p.lSp in 
powers of v. In particular, the (remarkable) fact that the (vanishing) IPN, 2PN 
(03) and 3PN (a^) contributions to the A potential are linear in v is already 
telling us something about both the "second-order GSF" (2 GSF) contribution 
i/^a2(u), and the "third-order GSF" (3 GSF) contribution j/^a3(w). Indeed, at 
the nPN order (which corresponds to a term oc in A{u)) one finds that, in 
intermediate calculations, the coefficient of is a polynomial in 1/ of degree 
n. It was, however, found in Ref. [25] that remarkable cancellations take place 
in the computation of the A function, and that the ly^ terms present at 2PN 
(oc u^), and both the i^^ terms and the terms present at 3PN (oc u^) exactly 
cancell in the final result for A^u;^). From a practical point of view, GSF 
studies have not yet embarked on any real "second-order GSF" work. We shall 
therefore focus, in the following, on the "first-order GSF" contributions v a{u) 
and iyd{u) in ([XTT)) . (|XT5)l . 

4 Circular orbits in the EOB formalism and the 
Last Stable (circular) Orbit (LSO) 

The stable circular orbits in the EOB formalism are conveniently discussed by 
using the dimensionless, rescaled variables (13. 8p and by considering the squared 
effective Hamiltonian for the radial motion, Eq. (|3.10p . As usual the presence of 
a positive kinetic energy term, / B+Q{pr), associated to the radial momentum, 
on the right-hand-side (RHS) of (|3.10p implies that the stable circular orbits (for 
a given dimensionless angular momentum j) correspond to minima (with respect 
to f) of the "effective radial potential" 

+ • (4-1) 

Using the short-hand notation (|3.14p . we shall then define the function 

W, [u; y) = A{u- y){l + ^ A{u- v) + fB{u; v) , (4.2) 

where we also introduced the short-hand 

B{u;v) =u^A{u]v) (4.3) 

which should not be confused with the metric component B[u] v) = entering 
(021). 

Stable circular orbits correspond to minima (with respect to u) of Wj (u), i.e. 
they solve W!j{u) — 0, with W'^'{u) > 0. [Here, and in the following, a prime will 
denote a u-derivative.] The solutions of Wj{u) — with Wj'(u) < correspond 
to unstable circular orbits, while the LSO is the solution of Wj{u) = which 
satisfies Wj'(u) = 0. 
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We can then parametrize the (one-parameter) sequence of circular orbits by 
the value oi u — 1/f. Indeed, while it is a priori difficult to solve (suppressing 
the presence oi u in A and B) 

W-{u) = A'{u)+fB'{u) ==0, (4.4) 

with respect to u, it is trivial to solve it with respect to j^, namely 

Knowing u and jcirc(w) one can then compute all the physical quantities 
attached to the circular orbit. From (|3.10p (and pr = 0) the (specific) effective 
energy is 

- ^A{u)+f^a^)Biu) , (4.6) 
while the corresponding total energy (divided by the total mass M ) is 



TTCITC j 

h{u) ^ ^ ^1 + 2K4T(^) - 1) . (4.7) 

Finally the orbital frequency around the circular orbit is obtained from 
Hamilton's equation of motion for ip: 

This yields 

n"'%u) = GM n'"^u) = ^"'■^(^)-^(") . (4.9) 

h{u) Hcs{u) 

Squaring (|4.9p . and inserting (|4.5p - (|4.7p in the result (remembering that 
B{u) = u^A{u)), finally leads to the following simple result 

»gi,,(^) _ lA'ju)^ 1 A'ju) 

2 h^u) 2 1 + 2z/(i/;:'f'-=(u) - 1) ■ 

The resuh (|4T0| is equivalent to Eq. (11) in [37 . Replacing tl = GM n by the 
corresponding x{Q) = {GM ^^/^ , Eq. (|2.18p . we can rewrite (|4.10p as 



^circ(^^) = U [ ^ ^ .ify. 1 , (4.11) 



1/3 

l + 2Ki?eT(«)-l) 

where the explicit expression of H"g'^{u) would be (from (|4.5p and (|4.6p ) 



rrrirr/ N /./ ^ ^'(^*) x / 2u A"^ (u) 
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Up to now we have made no approximations. In particular, we can conclude 
that the exact LSO frequency parameter xlso — (GMfiLSo)^^^ is obtained by 
inserting in (|4.1ip the root mlso of 

= A{u) EE A'{u) B"{u) - A"{u) B'{u) , (4.13) 

which is the condition for havine; a common solution to = A' + pB' = 
and W'j' = A" + pB" = 0. A more explicit form of the discriminant A is 

/\:^2AA' + Au{A'f -2uAA" . (4.14) 

As a simple check on the above results, and as a warm up for later, one easily 
sees that, in the test-mass limit — > 0, i.e. when A{u) ^ 1 — 2u, one recovers 
well-known results for a test mass in circular orbit around a Schwarzschild black 
holdZl 

3,Uu) = , , + 0{u) , (4.15) 

\/u(l — 3u) 

1 — 2?; 

H^J^u) = ^== + OH, (4.16) 
Vl — 3u 

{GM n,i,,iu)f/^ = Xcn-ciu) =u + 0[v) . (4.17) 
As for the LSO condition (|4.13p it reduces in this limit to 

= A(u) = 4(6u- 1) + 0(1/), (4.18) 

so that one recovers the well-known ulso = GM/rLso = ^ + ^('^)- 

One should also note that the full dynamics of circular orbits is entirely de- 
scribed by a single function, namely the EOB (— .gfjo ) uietric component A(u\ v). 
In particular the value of the LSO frequency parameter a;LSo(i^) can only depend 
on A(u\ v). 

Let us now insert in the general results (|4.5p - (|4.12p the "GSF expansion" of 
the EOB A potential, i.e. the expansion (j3.17p of the function A(u\ v) in powers 
of the symmetric mass ratio v. Keeping only the first-order correction in v ("1 
GSF approximation" ) , the GSF-expansion of Eq. (|4.5p yields 

Jcirc(^^)^ . ^ {\-\vd{u)~\v .f^""^ , +g(^^)) , (4.19) 

^Ju{\ - 3u) V 4 4 u(l - 3u) / 

where we introduced the short-hand notation 

h{u) = u^a{u), (4.20) 

so that B(u\ v) = u^(l — 2m) -|- vhiu) + OiiP'). Concerning the GSF-expansion 
of -ffeff, the zeroth-order (test-mass) result (j4.16p is sufficient for our present 



^We indeed recall that the EOB formalism is constructed so that the — > limit yields 
the dynamics of a test mass in a Schwarzschild spacetime. 
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purpose because HcS enters the frequency parameter x, Eq. (|4.1ip . multiplied 
by a factor v. The GSF-expansion of Eq. (|4.1ip yields 

..„(l-l.<.»-H.(-l=|=-l)+0(.')). (4.21) 

or, reciprocally, 

" = ^ (l + ^ ^ a'i.) + I . - l) + Oi.^)) . (4.22) 

To determine the location of the LSO we also need the GSF expansion of 
the discriminant A, Eq. (|4.13p . The expansion of Eq. (|4.14p yields 

- i A{u) = 1 - 6m + 1/ d{u) + 0{iy'^) , (4.23) 
where we introduced the short-hand 

a{u) = a(u ) + -(lOu - 1) a'{u) + - u(l - 2m) a"{u) . (4.24) 
Solving the LSO condition A(u) = then immediately leads to 

ui^S0 = l[l + vai^) + 0iiy^)] , (4.25) 

where 

«(B) = «(^) + J«'(i) + ^«"(^)- (4-26) 

Inserting the result (I4.24p in Eq. (|4.2ip finally yields the EOB prediction for 
the frequency parameter of the LSO to first order in 

:^LSO = ^[l + c^°^^^ + 0(i^')] (4.27) 



with 
where 

and 



EOB _ E , ~/l 



= c^+a(i), (4.28) 



1-/1 , (4.29) 



a(i)^a(i) + -a'(i) + -a"(i). (4.30) 

[Note that a(l/6) differs from a(l/6), Eq. (jiTM]) . in the coeflicient of a'.] The 
parameter c^°^ measures the fractional sensitivity of xlso to around = 0. 
It is also related to the fractional i/-sensitivity cff^ of the M adimensionalized 
LSO orbital frequency GM ^lso, 
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by 

= (4.32) 

As we see from (I4.28p . c^*^^ is made of two separate contributions: (i) a numer- 
ical contribution which is independent of the function a{u) (and which comes 
from the specific EOB "energy map" (|3.1[) relating the "effective" energy to 
the "real" energy); and (ii) a contribution which is a hnear combination of the 
values at the unperturbed LSO of a{u) and its first two derivatives. 

Let us now compare the EOB prediction (|4.28p to the recent GSF result 
(|2.20p . We shaU consider here the case where a{u) is given by a simple power- 
law PN expansion 

a{u) = a„ u" = 03 u"^ + 04 + u'^ + + aj + . . . (4.33) 

which is the type currently assumed in the NR-EOB comparisons. [Below, we 
shall argue that the coefficients 05, ae, . . . must include a logarithmic dependence 
on u.] Inserting the expansion (|4.33p in Eq. (|4.30p . we get (if we neglect any 
"logarithmic running" of 05, Og, • . •) 

^(i) = ^ a„ 03 -I- 04 -I- 05 + fle + 07 + . . . (4.34) 



where 

a„ = (2n2-n+l)^. (4.35) 
6" 

The currently known coefficients an in the PN expanqion of A{u-, v) [25] (see 
Eq. (|XT5l) ) are 

as = 2, 

94 41 TT^ 

04 = — — = 18.687903. (4.36) 

The contributions of the currently known terms in the PN expansion of A 
are, numerically, 

cf = 0.038127, (4.37) 

53 = ^ = 0.148148, (4.38) 
29 

04 = 04 = 0.418171, (4.39) 

1296 4 , \ ] 



so that they add up to 



03 + ^4 = 0.604446 . (4.40) 



Comparing this (partial) result to the GSF estimate (|2.20p . we see that the 
3PN approximation to a{u) explains about 72.5% of the GSF result. We expect 
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that the missing 27.5% will be contributed by the higher PN contributions to 
a{u), coming from 05, ag, a7, . . .. More precisely, by considering the difference 
between (|4.34l) and (|2.20p . we should expect (when neglecting any logarithmic 
running) that 

23 67 23 

■ «5 + 77^777; ae + — — 07 + . . . = 0.22975(40) . (4.41) 



3888 " 46656 " 69984 

Eq. (I4.4ip encapsulates the knowledge about higher-order (4PN, 5PN, 6PN,. . .) 
contributions to the function a{u) = [d A(u; v) / d i'],y=o contained in the recent 
GSF work [5^. We can roughly estimate the order of magnitude of the successive 
terms on the left-hand side (L.H.S.) of (|4.4ip by evaluating the (inverse) radius 
of convergence p of the series (j4.33p from the a4-coef5cient, say p ~ ^4 = = 
2.079171. As p4 is quite close to 2, it looks reasonable to expect that a„ will 
roughly grow as e„ 2" where = ± is a sign. This would mean that the higher- 
order contributions a„ to c^°^ decrease as e„(2n^ — n + 1)3~". For n = 6, 
this yields ±0.092 which is about 11% of the total a;™", Eq. while for 

n = 7 this yields ±0.042, which is about 5% of cjf". Though the decrease 
with n is rather slow, we can hope that the result (|4.41[) could give us, within 
a rough approximation, a constraint involving mainly 05 and og. Numerically 
this constraint would then read 

as + 0.242754 ag ~ 38.84(7) , (4.42) 

or equivalently 

4.11940 as + ag ~ 160.0(3). (4.43) 

The most accurate current EOB models use an A{u; v) potential which contains 
only two free parameters, denoted 05 and ag, and which has the following prop- 
erties: (i) A^^^[u, V, as, ag) is a certain Pade approximant, namely the ratio of 
two polynomials in u of the form, A^'^^{u\ v) = [1 ^ rii u)/{\ + di u + ^2 + 
^3 -\- d^u^ + ds u^): where the coefficients ni , di , . . . , ds £^re rational functions 
of V, and (ii) the expansion in powers of v of A[u] v) is of the form p.l7p with 



a^°^ (u) = 03 + a4 + as -u^ + ag , (4.44) 
af°^(u) = 67U^ + 68U^ + 69U^ + 6iou^" + 6iiu", (4.45) 

where the coefScients 67, ... , hn entering a^^^{u) are polynomials in a3, a4, as, 
ag. 

We see on Eq. (|4.13p that the i^-linearized a{u) function associated to these 
EOB models is of the general form (|4.33p written above, but with the restriction 
that a{u) is a polynomial in u which does not contain powers of u beyond 
. In addition, our notation has been chosen to be consistent in that the 
parameters as and ag entering (|4.44[) can be identified with the corresponding 
parameters in (|4.33|) . [In both cases too 03 and 04 are defined by Eq. (j4.36|l .] 
Therefore, when considering the current EOB models, the GSF-LSO constraint 
(|4.4ip does reduce to the linear constraint on 05 and ag written in (|4.42|) or 
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(|4.43p . In addition, the work of [THITB] has shown that the constraint of having 
a good agreement between NR data and EOB predictions selects a long and thin 
(banana-shaped) region in the as, ae plane. It was found by Damour and Nagar 
that this thin "good fit" region is, to a very good approximation, located around 
a line in the 05, og plane along which A'^qq {uq; 1; as, ag) is constant (and equal, 
say, to its value at (as, ag) = (—4, 24) which is one of the good- fit points, lying 
near the leftmost part of the good-fit region) . Here the prime denotes as above 
a u derivative; the first argument mq = (0.1)^/'^ = 0.215443 is a value of the u 
parameter which approximately corresponds to the EOB-predicted (adiabatic) 
LSO frequency for the equal-mass case namely MQ^^^{l' = 1/4) ~ 0.1; while 
the second argument in A'^q^ is the value of u corresponding to the equal- 
mass case {ly — 1/4), which is the case where very accurate NR data [36] were 
used to "tune" the considered EOB model. By numerically constructing the 
analytically defined central good-fit line just defined, one finds that it intersects 
the GSF-LSO straight fine (14.4211 into a unique point located at 

a^~-22.3, aE;'~+252. (4.46) 

These "intersection" (n) values for the parameters (asjOg) happen to lie 
roughly in the middle of the banana-shaped good-fit region found in [Ml [T6] 
(indeed, the latter region extends from (as,ag) = (0,-20) to (05,06) = (—36, 
-1-520)). This suggests that the information coming from the v 1 GSF-LSO 
study is able to break the degeneracy among as and ag left after fitting EOB to 
equal-mass (j/ — i) NR data. 

At this stage, the result (|4.46p has only an indicative value. It is a hint 
that the "real" values of the coefficients of the (4PN and 5PN) 1/ and 1/ 
contributions to the EOB potential A{u\ v) lie near the "intersection" values 
(|4.46p . To confirm (or refute) this hint one will need several types of further 
studies. First, it will be important to firm up the GSF result (|2.20p . [In the next 
Section, we shall indicate ways of doing so.] Second, one will need to explore in 
more detail the comparison between comparable-mass NR waveforms and EOB 
predictions to see: (i) whether the vicinity of the values (|4.46p does indeed leads 
to a better agreement; (ii) to what extent the possible "logarithmic running" 
of as and ag affects the estimates (I4.46P (see below) ; (iii) whether the influence 
of higher-PN coefficients 07, . . . is indeed small; and (iv) whether the current 
Pade-based definition of the EOB potential A{u\ v) is a sufficiently accurate 
representation of the A potential for all values of v in the interval < < |. 
Indeed, for the moment, the most accurate tests of the Pade-constructed A 
potential have been obtained for the equal-mass case, v = 1/4, for which the 
NR data are the most accurate [36]. Though the NR/EOB comparisons for 
several unequal- mass cases have also shown a good agreement, the numerical 
accuracy of these tests is smaller than the one reached in the equal-mass case. 
For all those reasons, it does not make sense, at this stage, to indicate "error 
bars" around the values (|4.46p . 



17 



5 Small-eccentricity orbits in the EOB formal- 
ism 



In the previous Section we considered exactly circular orbits in the EOB for- 
mahsm and showed that the GSF determination of the j/-derivative (at the 
test-mass hmit = 0) of the LSO orbital frequency was giving us access 
to a combination of radial derivatives (evaluated at the Schwarzschild LSO 
wSchw ^ g-^/^Schw ^ i^Q-^ of ^Yie i^-derivative a{u) = [d A{u;v)/dv]^^Q of 
the basic EOB radial function A{u; i') = — 5§o • Here, we shall consider (within 
the EOB formalism) slightly non circular orbits and show that a comparison 
of their gauge-invariant characteristics to GSF data can give us access to much 
more detailed information about the EOB formalism: in principle, one should be 
able to measure a certain gauge-invariant function p{x), where p{x) is a certain 
(x-dependent) linear combination of the z^-derivatives of the two basic EOB ra- 
dial functions, A{u; v) = — 5oo ^'^d B{u] v) = +g°^^ entering the effective metric 
(021). 

As we shall focus, in this Section, on the "small-eccentricity limit" where, 
in the effective Hamiltonian l|3.5p . the radial kinetic energy p^/B{r) (which 
is proportional to the square of the eccentricity e) is much smaller than the 
"azimuthal" kinetic energy J^/r^, we are allowed to neglect the higher-order 
contribution Q{r,pr) in (|3.5p . because it is (at least) quartic in pr, and therefore 
of higher-order {0{e'^)) in the eccentricity. [This result holds whatever be one's 
definition of the eccentricity e in this relativistic context.] Neglecting Q in 
the original EOB Hamilton-Jacobi equation (|3.3p leads to the usual Hamilton- 
Jacobi equation (/i-^ + g^^g p^ — 0) describing geodesic motion in the EOB 
effective metric p.2p . It is then convenient to describe this geodesic motion 
within a Lagrangian formalism, rather than the usual EOB Hamiltonian one. 
Using (Stueckelberg's) proper time formalism, we can then use as starting point 
the quadratic action 

5 = y'dreff/:(^a;^(reff),^^^) (5.1) 

where 

1 „ , \ , dx^ dx 1 1-1 , , 

Indeed, the critical points of the action (|5.ip . (|5.2p are geodesies of g°^. The 
value of the Lagrangian £ is a constant of motion which can be (after variation) 
constrained to have the value —^j?, corresponding to 

The constraint (j5.3p corresponds to using as (affine) parameter t^s along the 
geodesic the effective proper time divided by the (reduced) mass /Lt, 

dreff = ^ , (5.4) 
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so that the critical value of the action (|5.ip coincides with the usual "square- 
root" one — J fj,dscs- [For simplicity, we by-passed here the justification of 
the quadratic action (|5.2p based on the (Polyakov-like) use of an independent 
"einbein" degree of freedom along the worldline.] 

Inserting the effective EOB metric (I3.2p . one finds that the explicit form of 
the Lagrangian (|5.2p reads 



The Lagrangian (|5.5p admits the following conserved quantities 

feff = - .J,^, ^ A{r)-^, (5.6) 

d{dip/dTcg) dTcS ' 

£ = -Ai^ (5.8) 



Inserting (|5.6[) and (|5.7p into (|5.8p then leads to the following equation ruling 
the radial motion 

£2^ = A(r)B(r)(^A_y + A(r)(^M' + ^) • (5.9) 

It is easily checked that the radial dynamics defined by (|5.9p is equivalent to 
the small-eccentricity limit of the (squared) radial effective Hamiltonian (13. 6p . 
[The conserved effective angular momentum J having the same meaning in p.6p 
and (|5.9p . and the conserved effective energy £cs being numerically equal to 
the effective Hamiltonian HcS-] Passing to the dimensionless rescaled variables 
p.Sp . together with a corresponding (dimensionless) rescaled effective proper 
time f, such that 

Scff GM ^ 

Teff = = T, (5.10) 

^J. fi 

we find ^ 

- A{r) m (^g) + A{f) (l + fj) • (5.11) 

The rewriting of (|5.1ip in terms of the useful EOB radial variable u = 1/f = 
GM/r, yields the form 

?2 _A{u)B{u) f duV .2 2n 



geff = [,4 [-^) +Aiu)il+fu'). (5.12) 

The result (|5.12p exhibits how the circular-orbit effective potential Wj{u), Eq. 
(|4.2p . which depends only on the EOB potential A{u), gets modified by a radial- 
energy term which involves the product of the two functions A(u) and B{u), i.e. 
the combination D{u) introduced in Eq. (|3.12p . 
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Small-eccentricity orbits are solutions u = u{t) of the constraint (|5.12p of 
the form 

u{T)^uo + euiiT)+Oie^), (5.13) 

where uq is independent of f, and where e is a small parameter measuring the 
eccentricity of the orbit. Inserting (|5.13p into (|5.12p . one finds, at order e°, that 
uq must be an extremum of Wj [u) = A{u){\ + v?) = A{u) + B{u), i.e. that 
uo is related to j by 



We thereby recover that u{f) — uq corresponds to a circular orbit (see Eqs. (|4.4p . 
(|43)) above). 

Then at order (and modulo 0{e^) corrections), one finds that 



where W" = d'^Wj{u)/du^ (keeping j fixed). The constraint (|5.15p shows 
that iti(f) undergoes harmonic motion, i.e. ui(t) — acos{uJrT + ipo). Here 
(Dr denotes the frequency (measured in f time) of the small radial oscillations 
described by m = GM/r = uq + eui('r) + 0{e^). One sees on (|5.15p that is 
given by 

'^'M = ^;^W,"(uo), (5.16) 

where we used the notation (|3.12[) . More precisely, as long as uq is a minimum 
of the effective potential Wj{u), i.e. as long as Wj{uo) > 0, one finds that 
o)^ is positive, which indeed corresponds to (stable) harmonic radial oscillations 
Ui{t) (X cos{6jr t + ipo) . However, when Wo is a maximum oi Wj (u) (W^"(uo) <0 
so that tD^ < 0) the perturbations 5u = mi(t) of the circular orbit u = uq are 
unstable, and grow exponentially: ui(f) — aexpdtJrl t) -|-/3exp(— |d)r| "t), where 
Icj^l = \/—Cj1. The dividing hue is the LSO where W-(uq) = 0, i.e. a)^(wo) = 0, 
in agreement with the criterion (I4.13|) used above. 

We cannot directly relate the (EOB) f-time (squared) radial frequency CjI to 
the radial frequency studied in the recent GSF study [26] because the latter 
quantity refers to a different time variable r, namely the proper time of the 
Schwarzschild background (of mass m,2) used in GSF work. In order to trans- 
form the result ()5.16|) into a physical, gauge- invariant result, the simplest is to 
consider the dimensionless ratio of two frequencies: the radial frequency cD^j and 
the corresponding f-time azimuthal (or circular) frequency, say Cl = dip/dr. [We 
use a double hat on f2 to distinguish it from the dimensionless coordinate-time 
azimuthal frequency fi'^"'^ defined in Eq. (14. 9|) .] Using the angular momentum 
conservation law (|5.7p . we see that, along the circular orbit u — uq, the f-time 
azimuthal frequency is given by 

iliuo)^4o^j4. (5.17) 
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Finally, dividing (|5.16p by the square of (|5.17p leads to 
QjHuo) _ 1 1 fA"{uo) 



2 D{uo) ypiuo) 



B"iuo) , (5.18) 



where j^(wo) is given by Eq. (j5.14|) . Inserting the latter equation yields the 
explicit form 

D(uo) A{uo) 




(5.19) 



where A(u) is the determinant (|4.13p introduced above, and where D{u) de- 
notes, as in p.l3p . the inverse of D{u). 

The quantity on the L.H.S. of (|5.19p is gauge- invariant, and independent 
of the time parametrization. In particular, it is also equal to the square of 
the ratio of the coordinate-time radial frequency, say w^, to the coordinate- 
time azimuthal frequency = dip/dt. The R.H.S. is expressed in terms of the 
EOB (inverse) radial coordinate uq = GM/ro (along circular orbits). However, 
Eq. (j4.10p has related uq to the dimensionless, gauge-invariant frequency pa- 
rameter ricirc(uo) = GM fl'^"'^{uo) = GM d(p/dt. Therefore, the combination 
of (|4.10p and (|5.19p relates the two gauge-invariant observables (ujr/fl)'^ and 
GMQ. 

The EOB predictions (|4.10p and (|5.19p are valid for an arbitrary mass ratio 
u. Let us now focus on the limit v <ti 1, in which GSF studies can, in principl^, 
compute the L.H.S. of (|5.19p as a function oi ft ~ GM il. 

Inserting in (|5.19p the GSF expansions (|3.17p . p.lSp one finds (suppressing 

for simplicity the hats, as LOr/^ = cor/i^: the ratio of coordinate-time frequen- 
cies) 



V n J 



(1 - 6u) d{u) + a{u) + 2ua'{u) + ^u{l~ 2u) a" {u) 



= 1 — (!)U + V 

+ 0{v^), (5.20) 



where we henceforth also suppress the index on uq. By further inserting in 
(|5.20p the GSF expansion of u in terms of x, Eq. (|4.22p . we get 

[^)^ = l-&x + up{x)+0{v^), (5.21) 

where the function p{x) is made of three pieces: 

p{x) = pe{x) + pa{x) + pd{x) . (5.22) 



^Though there are some technical compHcations for deahng with unstable orbits (e.g. with 
extracting the conservative part of the dynamics), GSF studies can, in principle, compute the 
L.H.S. of | [5A9t for < -uo < 1/3, i.e. for ro > 3GM. 



21 



The first piece, pe{x), comes from tfie Energy map (|3.ip relating the EOB 
effective energy ScS — H^s to the EOB real energy Srcai — ^eob- It is given by 



Mx)=4.(^l--^j. (5.23) 

The second piece, Pa{x), is related to the function a{u). It reads 

p^[x) = a{x) + a; a'(x) + i a;(l - 2x) a" (x) . (5.24) 

Finally, the third piece, pd{x), is related to the fmiction d{u), and reads 

pd{x) ^ {I - 6x) d{x) . (5.25) 

Eqs. (I5.2ip - (I5.25|) represent one of the main new results of the present work. 
They show that GSF studies can give us access to a function of a; = G{mi + 
7712) which, when interpreted within the EOB formalism, is an x-dependent 
combination of a{x), a'{x), a"{x) and d(x), where a(x) and d{x) are the v- 
derivatives of the two basic functions A{u; v) and D(u; v) parametrizing the 
EOB effective metric. In principle, GSF calculations can numerically compute 
the function p{x) on the full interval < x < 1/3 corresponding to 3GM < 

To < +00. 

Note that the contribution Pd{x) contains a factor 1 — 6a; which vanishes 
at the (unperturbed) LSO. This implies that the value of the (z^-perturbed) 
LSO frequency obtained by requiring the vanishing of [uj^/^Y equivalent to 
that obtained by discarding the Pd{x) contribution to p{x). More precisely, the 
solution of {tOr/^Y = is 

a^LSO = \[l + vp{\) + 0{v^)\ 

= \[1 + HPE{\)+Pa{\))+0{V^)\ , (5.26) 

which is easily seen to coincide with our previous result (|4.27p - (l4.30p . 

On the other hand, though the exact (z^-linear) result ()5.2H) is equivalent 
to the exact (z^-linear) result ()4.27p - ()4.30|) for determining the exact (v-Ymeai) 
s^LSO J the PN expansion of (I5.2ip leads to a different estimate of the PN expan- 
sion of xlso than the PN expansion of (I4.27p . Indeed, let us consider the PN 
expansion (i.e. the expansion in powers of a;) of the exact function p{x): 

p^^ix) = p™(x) + p^ix) + pTi^) , (5.27) 

where 

pl^ix) =2x'-lx'-^x'~^-^x'-^-^x' + Oix') , (5.28) 

p^^{x) = 3a3x2 + (6a4 - 2a3)x^ + (lOas - 7a4)a;'* 

+ (15a6-14a5)a;^ + (21a7-23a6)a;'^ + C'(a;^), (5.29) 
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+ {d5-6d4)x^ + {de-6d5)x^ + 0{x'^). (5.30) 



Here the coefficients a„ (n > 3) are the expansion coefficients of the PN expan- 
sion (|4.33p of the function a{u). Similarly the coefficients dn {n > 2) are the 
expansion coefficients of the PN expansion of the function d(u), i.e. 

J(u) = ^J„u". (5.3f) 

n>2 

As above, we are here assuming, for simplicity, that the expansion coefficients 
a„ and dn are pure numbers (i.e. we neglect any possible logarithmic running 
starting at the 4PN level: 05,^4). Combining Eqs. ((O8)) - ((O0l) one finds for 
the PN expansion coefficients of p{x), i.e. 

p^^{x) = p2X^ + p3 x^ + p^x* + P5 x^ + P6 x^ + . . . (5.32) 



P2 = 


2 + 3a3 






P3 = 


2 




27 . 


P4 = 






675 


P5 = 


^"32" " 




3969 


P6 = 


64 



+ J2 , (5.33) 

13 + 604 — 6(12 + ^3, (5.34) 
' a4 + 10 as — 6 (is + c?4 , (5.35) 

14 05 + 1506-6^4 + 4, (5.36) 
-2306 + 2107-64 + ^6. (5.37) 

The 3PN-level calculation of the dynamics of two-body systems [38] has 
given us access to the values of 03 and 04 (see ()4.36p ). as well as of J2 and J3 
[29] . namely 

J2 = 6 , J3 = 52 . (5.38) 

Note that [25^ has also given us access to the beginning of the PN ex- 
pansions of the i/-quadratic (2GSF) (and even i/-cubic, 3GSF, for A{u; v)) 
terms in the GSF expansions ()3.17p . ()3.18|) of A{u;i') and D{u\v). Namely 
af^{u) = af^{u) = 0, and df^{u) = -6^^ (see Eqs. ({31^ . ^^). [We 
recall that the 3PN expansion of A{u) goes through the order u^, while that of 
D{u) stops at order ] 

Using the knowledge of 03, 04 (Eqs. (|4.36p ) and 3,2, d^ (Eqs. (|5.38p ) we can 
compute the first two terms in the PN expansion (j5.32p of p{x), namely 

P2 = 14, 

397 123 

= —^2 ^ 122.627416. (5.39) 

2 Id 

We have checked that the corresponding IGSF 3PN- accurate expansion of the 
ratio {uir/^Y , (I5.2ip . agrees with the (full GSF, 3PN-accurate) results derived 
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by Damour, Jaranowski and Schafer [35] ■ [Note that one must use the values 
Wkinotic = 41/24 [ini m] and = [38] in the results of [39].] Indeed these 
authors determined the PN expansions of both 

LOr = ^^radial = 27r/ P 

and 

^ ' — ^circ — ^radial ^pcriastron ^TT ^ 

The ratio (wr/^^)^ corresponds, in their notation, to (1 + fc)^^. For circular 
orbits, Ref. [39] computed the expansion in powers of 1/j^ of fc (and (1 + fc)^"*), 
see their Eqs (5.25)-(5.28). Inserting in these expansions their result (5.8) for 
the X expansion of j'circ one finds the i^-exact, 3PN-accurate expansion of (w^/f^)^ 
in powers of x. Its ^/-linear piece is found to agree with our result (|5.39p above. 

If one were to use the 3PN expansion of p(x) to estimate the LSO frequency, 
i.e. from the first equation (|5.26p (keeping the d-contribution), one would find 

a:LSo[p'™] - ^ [1 + V + 0{v^)\ (5.40) 

6 

with 

c.[p3PN]^g + g^ 0.956608. (5.41) 

If we compare this estimate to the previous 3PN-based estimate (|4.40p (based 
on the PN expansion of the exact, a(M)-dependent, result (|4.27p ). as well as to 
the recent GSF estimate (|2.20p . we see that (|5.4ip is now larger (by 14.67%) 
than the GSF estimate, and larger (by 58.26%) than the previous a?"^^ estimate 
(j4.40p . This shows the unreliability of using (non resummed) PN expansions for 
estimating physical quantities in the strong- field regime (here the LSO, which is 
a non-perturbative phenomenon). Although the p^^^ estimate (|5.4ip happens 
to be closer to the GSF result ^^W^ than the a^^^ one (jiliil)) . we think that 
this is purely accidental. Indeed, the results of Section 4 have shown that c^*^^ 
only depends on the a{u) function, and that the present 27.5% disagreement (on 
the low side) between C2:[a'^^^] and c^^^ is most probably due to neglecting the 
higher PN contributions to the a function, see Eq. (|4.4ip . By contrast, when 
comparing the various terms contribution to C2,[p'^^'^], ()5.41|) . one finds that it 
contains a rather large contribution connected with the large value of da which 
"pollutes" p3, see Eq. (|5.34|) . Indeed, note that the individual d-contributions 
to Ca;[p'^^^] are 



% ' (5-42) 



Pd2 7^ + Pd3 = S = S = 0.24074074 . (5.43) 





1 


d2 


Pd2 


62 


62 ' 




1 




Pdi 


63 


62 



so that their sum yields 
1 

62 ' 63 63 63 
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The crucial point is that while ^2 cancells (as it should) between the 2PN 
and 3PN contributions to Ca;[jO'^^'^], da = +52 does not cancell, and thereby 
artificially overestimates Cx- One sees on Eq. (|5.35p that the 4PN contribution to 
Cx \p^^^\ will introduce another contribution proportional to da which precisely 
cancells the contribution (|5.43p . Similarly one sees on Eqs. (|5.33p - (|5.37p that 
the contribution cx that appears at the 4PN level, Eq. (|5.35p . is cancelled 
by a contribution cx —6^4 entering the 5PN term (|5.36|) . In other words, while 
the exact expression (|5.25l) of pdix) shows that d[x) does not influence the LSO 
frequency, the PN expansion (|5.30p oi pdix) artificially introduces a contribution 
of d[x) at any finite PN order. 

A more general way of seeing the unreliability of using (non resummed) PN 
expansions to determine the LSO consists in considering, given any function of 
X (which does not vanish at the LSO), say /(x; v), the quantity 

F{x-v)^!{x-v){^\ (5.44) 

The exact solution of the constraint F(x) = is always the exact LSO, inde- 
pendently of the multiplicative factor /(x). However, if one considers the PN 
expansion of F{x) and solves some PN-truncated equation, say F'^^^{x) — 0, 
its solution will depend on the coefficients entering the PN expansion of ]{x). 

Going back to the issue of the knowledge that can be extracted from GSF 
calculations of the function [ujrl^'^ix)^ i.e. the knowledge of the function p(x), 
see Eq. (|5.2ip . we note that it would be very interesting to study this function 
not only near the LSO {x ^ ^), but also for small values of x (i.e. large radii). 
In this limit, one should be able to: (i) check the validity of the first two terms, 
P2 + P3 in the PN expansion of p(x\ and (ii) reliably extract some of the 
higher-order terms in the PN expansion ()5.32p (and explore their logarithmic 
running). Note in this respect that, in view of Eqs. (|5.33p - (|5.37p . measuring 
higher and higher coefficients p„ will give us access to some combinations of 
the higher PN coefficients of a{x) and d{x). More precisely, we will have the 
structure 

P4 10 05 -I- d4 , 

— 14 05 + 15 flg — 6 (i4 -I- J5 , 
P6 -23 06-1-2107-6 4-1-^6, (5.45) 

where the R.H.S.'s indicate the knowledge acquired from knowing the L.H.S.'s. 
Though it is somewhat frustrating to see the (less interesting) higher-order J„ 
coefficients entering (and "polluting") the R.H.S.'s we might still be able to 
reach some plausible results by combining (|5.45p with our previous result (|4.4ip 
which only contained the (more interesting) a„ coefficients. In particular, as it is 
easily checked that the main contribution to the 3PN coefficient p3, Eq. (|5.39p . 
comes from 04, it is plausible (especially in view of its rather large numerical 
coefficient) that the contribution 10 05 dominates in the quantity 10 05 -I- (i4 
which can be deduced from a measurement of p^. Then a GSF measurement of 
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P4 would give us an approximate knowledge of 05 . Inserting this knowledge in 
the approximate equation (|4.43p above, we would then be able to approximately 
determine uq too. It will then be interesting to see whether the values of 05 
and ae so determined are close to our EOB-NR-based estimate (|4.46p . If this 
is the case one might continue the approximation and try to deduce, e.g., an 
approximate estimate of 07. In addition, GSF data on the small x behaviour 
of p{x) will allow one to explore the presence of logarithmic terms starting at 
the 4PN level (see Section 7 below where one similarly discusses how GSF data 
on Detweiler's redshift function u^(x) can explore the presence of logarithmic 
terms in the latter function). 

6 Zero-binding zoom-whirl motion and other 
ways of extracting EOB information from 
GSF computations 

There are several other ways in which GSF computations might be used to 
extract information about the functions A{u), D{u), Q{u,pr) entering the EOB 
formalism. Let us here sketch a few possibilities. 

First, one might go beyond the small-eccentricity limit and compare IGSF 
computations of eccentric orbits to EOB predictions. However, to do so one 
should include the effect of the 0{pf) contribution Q{u,pr)- At present one 
only knows the 3PN-level expression (|3.11|1 of Q{u,pr) ^ u^p*. At higher 
orders it will probably involve both higher powers of u and higher powers of Pr , 
say 

Qiu,pr) = ,^[q{u)pt + qiu)pt] + 0{v^) + 0{pl) , (6.1) 

with q{u) = (72 + 93 u'^ + . . .. A detailed EOB/GSF comparison for smallish- 
eccentricity orbits might allow one to probe the function q{ii), together with 
a{u) and d{u). 

Second, a more ambitious project might be to try to extract the total con- 
served (EOB) energy and angular momentum, £ and J from GSF computations. 
In principle, one GSF way of determining £ and to the required accuracy 
would be to go to the second GSF approximation (2GSF), i.e. to compute the 
metric perturbation through order v"^ , and then to compute, from the metric, 
the ADM surface integrals giving the total energy and angular momentum of the 
system. This will give £-Af = ei v + e2 iy'^ + 0{iy^) and J = ji i'+j2 i''^ + 0{iy^), 
which is the precision needed to explore the effects linked to i^a^u), vd{u)^ . . . 
However, it is not clear when the GSF community will be in position to compute 
the 2GSF approximation. Let us therefore sketch what can be done now with 
the IGSF approximation. The knowledge of the IGSF conservative self-force 
allows one, in principle, to determine some conserved energy-like and angular- 
momentum-like quantities along perturbed orbits. The clearest way to do so is 
to consider unbound orbits, that come from some "in state" with infinite sepa- 
ration between mi and TO2. The "in state" then determines, in principle, the 
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numerical value of £ and J. One must, however, be careful to include the effect 
of the recoil of the large mass m2- In an Hamiltonian formalism (such as the 
ADM or the EOB one), the total conserved angular momentum is given by 



(6.2) 



while the total conserved energy has the form 



£ = H 



El 



(6.3) 



where E'^ = {rn^ +p^)^/^, a — 1,2, denotes the free kinetic mass-energy, and 
where the interaction term Hint tends to zero in the infinite separation limit 
r — > 00. In addition, the conditions to be in the center of mass frame have the 
form 

P1+P2-0, (6.4) 
E'I Xi + E2 X2 + interaction terms — , 



(6.5) 



where the interaction terms are proportional to the gravitational constant, and 
are smaller by a factor r than the leading terms in (|6.5p (see [lOj)- Usually, 
GSF calculations only follow the dynamics of the small mass mi. It would then 
record (for unbound orbits) only the mi contribution to £ and coming from 
the incoming state (with infinite separation), say 



= [Xi X pi]c 



(6.6) 



mi 



'Pi 



(6.7) 



In view of (|6.2p - (|6.5p one must then "renormalize" both the angular mo- 
mentum and the energy 



f^oc 

J = 1 1 + =^ ) jr 



£ = Ef 



Pi 



(6.8) 



(6.9) 



Having so computed £ and J one might then, for instance, compare the 
EOB prediction for the scattering angle 9{£^ J) (which follows from the EOB 
Hamiltonian) with GSF computations of 6 for a sample of values of £ and J'. We 
see that, in principle, we have here access to one function of two real variables, 
which is ample information for determining the functions entering the EOB 
formalism. 

We have been assuming here that, with some effort (notably concerning the 
separation of the conservative part of the self- force from its radiative part), one 
could deal with unbound orbits. In case it is computationally much easier to 
deal with bound orbits, one might still be able to extract £ and J from GSF 
computations. First, we remark that in view of the parity properties of the 
conservative self- force under time reversal (see, e.g., [26|) the ratios Rt = Ft/u^ 
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and Rip = Fi^ju^ (where vJ' = dr/dr, using the GSF proper time of, say, [55]) 
are even functions of r. [Here, we fix t = at an extremum of the radial 
coordinate.] This even-parity property is shared by the radius as a function of 
T. We can then follow a strategy used in [55] and (numerically) consider that 
the ratios Rt and R^ are some functions of r. 

This leads to evolution equations for the "toi contribution" to the energy 
and angular momentum of the form 

±{m,Ut)^F,^^R,{r-ElLl), (6.10) 

^{miu^) = F^ = ^R^{r-El,Ll). (6.11) 
ctr ctr 

Multiplying both sides of these equations by dr one sees that the consideration 
of the radial integrals / dr Rt (r) and / dr R^ (r) allows one to "correct" the 
"bare particle energy" —mi ut and the "bare particle angular momentum" by 
additional functions of r so as to define some IGSF-conserved "improved" par- 
ticle energy and angular momentum. By exploring the r-dependence of these 
improved quantities as r gets large, one might be able to so deduce the quan- 
tities and , Eqs. (|6.6[) . (|6.7p . i.e. the values of the mi contributions for 
infinite separation. It then remains to "correct" them so that, as in Eqs. (|6.8[) . 
(|6.9|) above, they also incorporate the "recoil" contribution (or, more precisely, 
the analytic continuation in £ and of the functions that express these recoil 
properties in the infinite-separation limit). Finally, if one succeeds in determin- 
ing £ and J for bound orbits, the consideration, say as in Ref. [31], of the two 
gauge-invariant functions of two gauge-invariant variables 

^radial = '^r(f , J) , (6-12) 
"^pcriastron — 

{£,J), (6.13) 

will give us access to ample information for determining the functions a(u), diu), 
Q{u,Pr) entering the EOB formalism. 

Let us end this Section by considering a special motion which should be 
easier to investigate by GSF means, and which should give us access to very 
useful information about the crucial EOB function a(u), independently from the 
other ingredients of the EOB formalism. Specifically, we wish to consider here 
the special zero-binding zoom-whirl orbit which starts, in the infinite part, with 
zero kinetic energy (but a non-zero angular momentum) at infinite separation 
and ends up, in the infinite future, "whirling" indefinitely around some limiting 
finite separation. In the Schwarzschild case {v limit) this motion has 
dimensionless angular momentum = 4 and its limiting dimensionless whirl 
radius is fwhiri = 4. Our purpose is to study the EOB predictions for the 
modifications of these values when ly ^ 0. From the GSF point of view, once 
one has a way to compute the conservative part of the self-force (for such an 
unbound orbit), it should be straightforward to determine which deviation from 
Ui^ = 4 is needed to end up whirling undefinitely around the large mass. Having 
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so determined both (taking into account the "correction" (|6.8p l. and f whirl, 
one will also have access to the final orbital frequency f^whiri corresponding 
to fwhiii- In other words, this special zero-binding zoom- whirl motion gives 
us access to two dimensionless observables: and fi^hiri- Then, as discussed 
below, one can also measure a t/iird dimensionless observable: the value oiuj^/Vi^ 
at the whirl radius. 

From the EOB point of view the special motion we are considering has a 
total energy £ — H^ob equal to the total mass M. From Eq. (|3.1[) we see that 
this "zero-binding condition" is equivalent to requiring that the dimensionless 
effective energy £off — H^s be equal to 1 . In the EOB formalism a generic radial 
motion is described by the constraint 

= 4ff = const, (6.14) 
where u= 1/f, and where 

Wj{u) = A{u){l+fu'^) (6.15) 

is the "effective radial potential" . Among the motions described by (|6.14p we are 
interested in the special case where f starts, in the infinite past, at f = -t-oo with 
zero radial velocity (and therefore zero radial momentum pr), and ends, in the 
infinite future, at some limiting radius rwhirb with zero radial velocity again. It 
is easily seen that these two conditions amounts to saying that the "horizontal" 
level line = 1, in the plane (f,f^g) [or {u,£es)], must satisfy the following 
conditions: (i) it intersects the "effective potential curve" = Wj{u), both 
at f = oo and at f = Twhiri; and (ii) it is tangent to the effective potential 
curve at f = f^hiri- The condition (i) expresses the fact that f — oo (or u = 0) 
and f = fwhiri (or u = tiwhiri = l/?^wiiiri) are turning points of the radial motion 
{Pr = 0), while the condition (ii) expresses the fact that the formal turning point 
f = ^whiri is reached only after an infinite time. In other words, (ii) expresses the 
fact that f ~ ^whiri is an unstable circular orbit. Note also that, from Eq. (j3.ip . 
the condition £cff — 1 does correspond to £ = i?EOB = M { "zero binding" ) . In 
terms of equations, (i) and (ii) imply that the special value j» corresponding to 
the zero-binding zoom-whirl motion must satisfy 

Wj, (Uwhirl) = ^(Wwhirl) + J* ^(Wwhirl) = 1 , (6.16) 
(Wwhirl) = ^'(Wwhirl) + £ S'(Uwhirl) = . (6.17) 

Note the fact that the conditions (|6.16p . (|6.17p for the special zero-binding 
zoom-whirl motion depends only on the A(u; v) function (remembering that 
B{u) is just a short-hand notation for u'^A{u)). In particular, it does not depend 
on the higher-order Q{u,pr) term in (|6.14p . Indeed, though we do not know 
this term beyond the 3PN approximation (|3.1ip . it is enough to know that 
Q(u,pr) exactly vanishes as 0(j)f) as Pr —>■ 0. [Actually even a weaker vanishing 
oc p^ would suffice.] This remarkable property of this special motion makes it 
especially worth to be explored by a GSF approach. 



W,{u)+A{u) 



Bin) 



Q{u,pr) 
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Eliminating from these two equations (which is equivalent to using the 
circularity condition (|4.5p ) leads to the condition 

AB'-A'B = B'. (6.18) 

Remembering that B{u) is just a short-hand notation for A{u), the condition 
(j6.18p can be written down more explicitly as 

A^=A+^uA'. (6.19) 

In Eqs. (|6.18p and (|6.19p we have provisionally suppressed the mention that the 
argument u must take the special value u = Uwhiri- 

So far we have made no approximations. Let us now consider the IGSF 
approximation to the exact conditions (|6.17p . (|6.19p . i.e. let us insert in them 
A{u; v) = l-2u + y a{u) + Oiv"^). The IGSF expansion of ((6l9l) reads 

m(4m- 1) + (1 - 4u)i/a(u) = ]^vua'{u). (6.20) 

In the V ^ Q limit one recovers that its solution is Uwhiri = i + ^{^): 
mentioned above. Inserting this zero-th order knowledge back in (I6.20p . one 
finds that the IGSF (i^-linear) modification of MwWri is predicted, by the EOB 
formalism, to be 



11 

^whiil — ~ T 
''"whirl ^ 



l + -a'{\)+0{u')\ . (6.21) 



By itself, this result is not gauge-invariant as it refers to an EOB radial co- 
ordinate. Let us now convert the result (|6.2ip into physical, gauge-invariant 
results. 

On the one hand, by inserting (|6.2ip into either (|6.16p or (|6.17p one can 
compute the special total angular momentum of the zero-binding zoom- whirl 
motion. One finds 

j,^A[l~2iyail)+Oii^^)] . (6.22) 

In addition, by inserting (j6.2ip into our previous result (|4.2ip . we can com- 
pute the ("flat coordinate time") orbital frequency parameter of the limiting 
zero-binding circular orbit 



_ 1 

^whirl 



l + ^a'(i)+0(^2)] . (6.23) 

3/2 

The corresponding dimcnsionlcss frequency fiwhiri — -^^^whiri = ^^whiri 
then given by 

i^whiri = Ml + ^ + Oiiy^)] . (6.24) 

Note the simplifying facts that j*, Eq. (|6.22p . depends only on the value of 
the a function at u = 1/4, while a^whiri and ilwhiri, Eqs. (|6.23p . (|6.24p . only de- 
pend on the value there of the derivative of a(u). [The "energy terms" involving 
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(1 — 2u)l\J\ — 3u— 1 in (I4.2ip have also disappeared because of the zero-binding 
condition.] Finally, by extending the calculations of [26] to unstable circular or- 
bits, one might also have access to the value of the ratio w^/fi^ (with < 0) 
at our special zero-binding limiting circular orbit (|6.2ip . In other words, one 
might also, using Eq. ([OT|) . know the value of Using Eqs. 

and our results (|6.22p , (|6.23p , this will give us the value of the following combi- 
nation, — \ involving now the second derivative of a(u) and the 
function d[u). Therefore a GSF computation of the characteristics of the special 
zero-binding zoom-whirl motion would give use access to two interesting pieces 
of information concerning the strong-field behaviour of the ai^i) function: a(-|) 
and a'(-j). In addition, it will give us the value of a combination of a"(|;) and 
As the argument ~ 1/4 is larger than the (unperturbed) LSO value 
ulso — 1/6, the knowledge of a(-|) and a'(^) will provide a strong constraint 
on the shape of the function aiu) for all the arguments u < \. Note also that 
the PN expansion of a(-|) reads (neglecting any logarithmic running) 

iW 0,3 a-j 
°(4) = ^ + ^ + ^ + ^ + ^ + --- (6-25) 

If the rough estimate mentionned above, a„ ~ ± 2", holds this series would 
behave as ^ ± 2~". In particular, the n — 1 contribution might be as small 

n 

as ±1/128, i.e. smaller than 1%. Then we can use the current 3PN results 
(|4.36p . together with our tentative estimate (|4.46p of 05 and ag, to "predict" an 
approximate estimate of a(-|), namely 

a{\) ~ 0.03125 + 0.07300- 0.02178 + 0.06152 

~ 0.144. (6.26) 

The first line of ()6.26|) indicates the successive contributions of 03, 04, 05 and oe. 
Because of the large value of Og , (|4.46p . which is somewhat larger than the naive 
expectation ±2^ = ±64 (indeed {a^fl^ ~ 2.51 is larger than a^'^ ~ 2.079), 
the series (j6.26l) does not exhibit a clear convergence. However, even if a„ 
grows more like (2.51)", one still expects the n — 1 contribution to be only 
±(2.51/4)7 = ±0.038. In other words, one might hope that a(i) = 0.144 ± 
0.040, corresponding to = 4 [1 - 0.288(80) z^ + ©(i/^)]. It will be interesting to 
see how this approximate prediction compares to GSF estimates. Note finally 
that the corresponding PN expansion for a'(i) has a worse convergence because 
(neglecting any logarithmic running) 

n>3 

includes a growing coefficient n. Because of the worsened convergence of (|6.27p 
it would be poorly justified to venture a numerical estimate of o'(|;), beyond 
the fact that one expects it to be positive (on the basis of 03 + 04 giving the 
contribution +1.54 and og + 07 giving +1.10). 
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7 About logarithmic terms in PN expansions 



In the previous Sections we have considered, for simphcity, that the PN expan- 
sions of the various EOB functions (such as A{u, u) or its v-Ymeai contribution 
a(u)) contained only powers of the PN expansion parameter (u = GM/c^r or 
X = (GM Q,/ c^Y^^ for functions of x). This simphfying assumption aUowed us 
to connect our results to the recent results of the NR/EOB comparison which 
have all assumed that the function A{u; v) admitted such power-law PN ex- 
pansions, bacause they define A{u) by Pade-resumming the truncated Taylor 
expansion: 

^Taylor = 1 - 2m -f v{aj, + + + u'^ + . . .) + 0{v'^) . (7.1) 

We wish, however, to point out that the PN expansion of A{u), contains log- 
arithmic terms, starting at the 4PN level. This means that the coefficients 
asjOg,... in Eq. (|7.ip are not numerical constants but exhibit a "logarithmic 
running" with u, probably of the type 

05 (In u) — flg -I- flg In w , 

a6(ln u) — ftg -I- a\ In u . (7.2) 

Indeed, we shall argue below that the logarithmic running of 05 and ag is linear 
in In u. We, however, expect that some of the higher-order PN coefficients will 
exhibit some nonlinear runnings, say 

a„(ln u) = a° + a,\ In it + ... + aP(ln uf , (7.3) 

where the maximum power p increases (not faster than n) with n (probably 
p<n - A). 

Actually, as the A{u) function is only related in a rather indirect (and "non 
local" ) way to the spacetime metric of a binary system, we shall discuss here 
the issue of logarithmic terms in the context of another function, defined along 
the sequence of circular orbits, which has been recently studied in detail within 
the GSF approach. We mean here Detweiler's gauge-invariant circular-orbit 
"redshift" function u*-{jj) [57], where y = (Gm2 fi/c'^)^/^ (which differs from the 
"symmetric" frequency parameter x by the replacement M = mi + m2 —f TO2)- 
Expanding u*'{y) to first-order in the mass ratio q = mi/m2 ^ 1 (so that 

u\y) - (1 - 3z/)-l/2 -qu{y) + 0{q^) , (7.4) 

defines a IGSF-level function u{y) (here defined with a minus sign, so that u{y) 
be positive and involve mainly positive coefiicients). Recently, Ref. [28J has 
determined the 3PN-accurate post-Newtonian expansion of u{y) (= —w^p{y) in 
their notation), namely 

u{y) =uiy + U2y^ + U3y^ + Th y^ + 0{y'^) , (7.5) 
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where ^ [55] 



ui = l; U2 = 2; U3 = 5 ; "4 = 32 ' ^'^'^^ 

Let us first point out that there seems to be a simple connection between the 
PN expansion coefficients u„ of u(jj), and the PN expansion coefficients a„ of 
the EOB function a(u). Indeed, we can write 

Un = a-n + k„ , (7.7) 

where oi = = 02 and the higher ones are given in Eq. (|4.33p . and where 

ki^l; fc2 = 2 ; fca = 3 ; ki ^ 9 . (7.8) 

It is remarkable that the difference u„ — a„ is an integer. [This might, however, 
be true only for the lowest PN orders, as A:„ might also depend on the lower- 
order ttm's (with TO < n).] Postponing to future work a detailed discussion 
of the link between u{y) and a{u), we shall content ourselves here to note the 
simple connection (|7.7p . (|7.8p and to heuristically argue that the presence of 
logarithmic terms in the PN expansion of u{y) will entail the presence of related 
logarithmic terms in the PN expansion of a{u). 

We now focus on the presence of logarithmic terms in the PN expansion of 
u{y), which is more directly related to the PN expansion of the spacetime metric 

of a binary system. The issue of logarithmic terms in PN expansions of g^i, 
has been studied some time ago [42l |43j |44] . Let us summarize here the main 
results. First, one should distinguish the issue of PN logarithms arising in the 
"external metric" (and also the "wave zone" metric), from that of PN logarithms 
arising in the "inner metric" , or "near zone" metric. Logarithms enter the PN 
expansion of the external metric at the 3PN level [421144,- They are linked to the 
cubic interaction (mass monopole) x (mass monopole) x (mass multipole) (see 
Eq. (A6) of pij). These terms enter the relation between the radiative multipole 
moments and the multipole moments of the source (and show up in the energy 
flux of binary systems at the 3PN level, see e.g. Eq. (231) in the review [4]). 
Here, we are interested in the logarithms arising in the PN expansion of the 
inner (near-zone) metric. It was shown in [44 that the first logarithmic terms 
in the PN expansion of the inner metric arise at the 4PN level, and read (in a 
suitable gauge@), from Eq. (6.39) there, 

where I{t) is the mass monopole of the source (which can be approximated here 
as I{t) ~ M = TOi -|- TO2), and ij^J the sixth time derivative of the quadrupole 

^Though some gauge transformations can be problematic (as discussed above for the Lorenz 
gauge), we have checked that the gauge transformations used in .44 do not affect the present 
computation. 
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moment of the source (which can be approximated here by its Newtonian esti- 
mate lab — iTT-i Ui^^^ + m2 y'^^^ , where , ^2 ^-re the positions of the two masses, 
and where the angular brackets denote a symmetric trace-free (STF) projection). 
Remarkably, though the logarithmic term (|7.9p comes from a (hereditary) "tail" 
correction to the leading (Burke-Thorne) radiation-reaction contribution to the 

inner metric (5oo'*~™'*'^ = — | /^^■'/c'', see [H]), the logarithmic term (|7.9p 
is symmetric under time reversal, and thereby survives as is in the conserva- 
tive dynamics of the system. [This is due to the fact that the full hereditary 
radiation-reaction term that gives rise to (|7.9p . is given by a "logarithmic tail 
integral" over the entire past (see Eq. (6.38) of |44j ) which is time- asymmetric 
without being time-odd as the leading radiation-reaction term.] 

Let us now consider the contribution of the logarithmic term ()7.9p to De- 
tweiler's redshift function (along the sequence of circular orbits of the "small 
mass" mi) 

C- \ -1/2 

where t/^^, = g^i^y{yi) denotes the (regularized) value of the spacetime metric 
at the location y\ of the "small mass" mi , and where v\ denotes its coordinate 
velocity dy\/dt. 

We note first that the logarithmic term (|7.9p is easily seen to be proportional 
to the symmetric mass ratio v, and therefore to contribute directly to the u{y) 
piece within u*{y). There are, a priori, two sorts of contributions to u*(j/) 
entailed by (17. 9p . a direct one 

Jdi-ct^t^^)^ l^^i„(y^) (7.11) 

coming from the —goo contribution in (|7.10p . and an indirect one coming from 
the fact that u*(y) should be computed for a fixed value of the orbital frequency 
fl, which implies that the coordinate radius ri of the orbit of y\ must be per- 
turbed by the additional term (j7.9p in the metric with respect to the value it 
had before taking this term into account. However, a simplification occurs in 
that the coefficient of the perturbation 6ri of ?'i (approximately) vanishes. In- 
deed, the "Newtonian" contribution to the expansion of the R.H.S. of (|7.10p is 
{ri2 = |yi2|, 2/12 =yi- 2/2) 

^ + U = ^ + l{^yn^rl,. (7.12) 
ri2 2 ri2 2 V Al / 

In the limit mi <C m2, the ri2 (or ri ~ ri2) derivative of (|7.12p vanishes by virtue 
of the Newtonian equations of motion, i.e. Kepler's law fi^ ri2 = GM. [This also 
follows from the fact that we are expanding here the inverse of the Lagrangian 
of the particle yi] This simplification means that Eq. (|7.1ip represents (in the 
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limit mi <^ 7712) the full contribution to —qu{y) in Eq. (|7.4p . In other words, 
we have 

q S''^' u{y) = +i^^y1y>lM^ y[f ) . (7. 13) 

Evaluating the time derivative in (|7.13p for circular orbits (in the center of 

mass of the system, so that = ^2/12''^ ^^'^ 2/i = {'m-2/M) yi2, as already used 
above), and noticing that y ~ {Gm2 ^/c'Y^'^ is proportional to \/c? so that In y 
contains —2 In c, we find from (|7.13p that the leading logarithmic contribution 
to u{y) must be 

fi4 

<5'°«jl(y) =+-2/'Mny. (7.14) 
5 

In other words, this is saying that the term formally denoted 0{y^) in 
Eq. (frsjl is of the form 

{ul + ul\ny)y^ (7.15) 

with u\ = 64/5 = 12.8. Using the numerical data reported in Table I of [28], 
we have confirmed that this seems to be indeed the case with ~ 12.8 and 
~ 114.4. 

By extending the arguments of [44] leading to (|7.9p . it seems that the next 
logarithmic terms will be ©(In c/c^^), i.e. at the 5PN level, and involving only 
the first power of In c instead of a possible (In c)^ (which formally enters higher 
tail terms). This suggests that the next term is the PN expansion of u{y) will 
be of the form 

iu°+ullny)y\ (7.16) 

where m| is an analytically calculable quantity. A fit of the numerical data of 
[28] suggests that is comparable to ul (i.e. ~ 12) and that Ug ~ 360. At 
higher PN levels, one expects higher powers of In y to arise. The growth of the 
power of In c with the iteration order was shown very generally in [43 and j44j 
(see Eq. (5.2) there). See also Eq. (5.3) in [44] which says that the power of In c 
(say p) grows linearly with the PN order (say N): p = N — c. The constant 
c was 3 in [4^, if we gauge the PN order by the power of 1/c in the spatial 
metric. The GSF data behind the figures given in [^8 should be able to explore 
these phenomena. From the point of view of the present paper, which explores 
the contacts between the EOB formalism and GSF results, our main (tentative) 
conclusion (based part on the "experimental link" (|7.7p and part on the fact that 
the A{u) function can be, in principle, computed from the spacetime metric) 
is the corresponding occurence of logarithmic terms in the EOB formalism, as 
sketched in Eqs. (|7.2p . (|7.3p above. We leave to future work a computation of 
the logarithmic coefficient 05 in (|7.2p . 

As a final comment let us mention a possible subtlety. The logarithmic 
terms we have discussed above are infrared (IR) logs (linked to the matching 
between the near-zone and the wave-zone). There could also exist additional 
ultraviolet (UV) logs, linked to "finite size effects" , i.e. to the matching between 
the "body zones" (of order the radii of the considered compact objetcs) and 
the near-zone. From an effective field theory point of view, logarithmic terms 
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linked to finite-size effects would correspond to (logarithmic) UV divergencies in 
a pure point-mass description of a two-body system. Some time ago [45] , it was 
shown that the effective description of two gravitationaly compact bodies as two 
point masses (with action Spointmass = —^A J mAdsA, together with analytic 
or dimensional regularization) is valid up to the 5PN level, where finite-size 
effects (linked to tidal interactions, and depending on the internal structure 
of the compact bodies via a certain "relativistic Love number A:" ; see p. 83 
of [l5]) start appearing. The finite-size effects of body 1 correspond to the 
appearance of a quadrupole term in the metric, given (at the leading, Newtonian, 
approximation) by (see Eq. (19) in Section 5 of [15] ) 

q6 ^ 
(^ffoo ^ -TT h ml 1712 di-j — di-j r . (7-17) 

[The label 1 on fc refers to body 1. We suppress here the unrelated index 2 on the 
Love number k referring to its quadrupolar nature.] The additional finite-size 
effect (|7.17p can be described by augmenting the action of the two point masses 
by a non-minimal worldline coupling of the form T,a j I^a J dsA £^0 , where 
£^ = [u'^ u'^ C^avp]^ is the "electric" worldline projection of the Weyl tensor 
[lolHT]. From an effective field theory point of view, one expects that the addi- 
tion of such non-minimal couplings is needed not only to describe the finite part 
of extension effects, but also to "renormalize" the divergences (1/e poles in di- 
mensional regularization) that arise when one tries to describe extended objects 
by a point-mass action. It is therefore a priori plausible that a dimensional- 
regularization (DR) calculation of the interaction of two point masses generate 
1/e poles at the 5PN level, linked to terms of the type (|7.17p . but with a co- 
efficient oc 1/e. As is well-known, 1/e poles in DR correspond to logarithmic 
divergences, and always come accompanied by some logarithm of the ratio of 
the two relevant length scales: here that of the body zone, and the scale of 
variation of the metric near the considered body. This argument therefore sug- 
gests that the logarithm of (body size)/(scale of variation) ~ GM j (? ryi-, i-S- 
In u can arise at the 5PN level as an IJV log. In other words, at the 5PN level 
{0{u^) in 500) there might be two sources of In u: IR and UV. On the other 
hand, the recent work on tidal effects in neutron stars [47l [48] has shown that 
the (quadrupolar) tidal coupling coefhcient fiA oc kA formally tended towards 
zero as the compactness of the neutron star tended towards that of a black 
hole. This may mean that the a priori expected 5PN finite-size divergencies of 
gravitationally interacting point masses cancell out, and do not give rise to 1/e 
poles. In this case, there will be no UV source of In u at the 5PN level, and the 
generalization of the 4PN IR argument above should give the entire 5PN logs. 



8 Conclusions 

We have discussed various ways in which the computation of conservative Grav- 
itational Self Force (GSF) effects on a point mass moving in a Schwarzschild 
background can inform us about some of the basic functions, A{u;i'), D[u]v), 
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Q{u,pr), entering the Effective One Body (EOB) formalism. [Here u = G{mi + 
m2)/r and v = mim2 / {mi + 7712)^.] Our main conclusions are the following: 

1. The recently published GSF calculation ^26j of the 0{iy^) shift of the 
orbital frequency of the Last Stable (circular) Orbit (LSO) gives us access to 
the combination (|4.30p of the values of a{u) = [d A{u;i')/di'],^^o, Eq. p.l7p . 
and its first two derivatives, a'{u) and a"{u), taken at the unperturbed LSO 
location ulso = ^ + C'(j^). 

2. The current 3PN-based [221 [3H] knowledge of the beginning of the Taylor 
expansion of the function a{u) "explains" 72.5% of the GSF result. We expect 
that the missing 27.5% will be contributed by the higher post-Newtonian (PN) 
contributions to the function a{u), Eq. (|4.33|) . and notably by (4PN con- 
tribution) and ae (5PN contribution). Combining the GSF result [26] and the 
3PN-E0B results [29l|38], we determined in Eq (|4.41|) the value of a hnear com- 
bination of a5, ae, 07, etc. . . If the contributions of 07, ag, etc. . . are relatively 
small (which can be argued for) this leads to determining a linear combination 
of as and ag, see Eqs. (|4.42l) or (|4.43p . 

3. By combining the just mentionned approximate knowledge of a linear 
combination of 05 and ag with the recent determination [TU [TB] of the long and 
thin region of the (05,05) plane where the most accurate current EOB models 
[ini m] exhibit an excellent agreement with NR data for comparable masses 
(4i/ = C(l)), we found that they intersected near the point (|4.46p of the (05,05) 
plane. This suggests that the information coming from the v 1 GSF study 
of the LSO is able to break the degeneracy among 05 and oe left after tuning 
the two-parameter EOB (05, 05) waveform to equal-mass {i^ = j) NR data [55] . 
More work (both on the EOB side and on the GSF one) is however needed to 
confirm the tentative values (|4.46p . [See the end of Section 4 for a more detailed 
discussion.] 

Note that the values (|4.46p . when inserted in the EOB formalism, predict 
the value of the LSO orbital frequency for all values of ly, between the test-mass 
limit 1^ <^ 1 and the equal-mass case ly = 0.25. The exact EOB prediction is 
obtained by using the exact EOB results given in Section 4 (notably Eq. (|4.9p V 
Let us only indicate here that the final result for the j/-dependence of the GSF 
n NR-tuned f^LSO ('^) can be approximately fitted by a quadratic polynomial in 

V of the form 

CIlsoH = G{mi + m2) nLsoH-G'^^^^il + 1-25 iy+ 1.87 ly^]. (8-1) 

It should be noted on the result (|8.ip that the non-linear dependence on 

V (here summarized by the term +1.87v^, but given in reality by a more 
complicated function 021^^ + 031^^ + ...) is numerically quite important for 
comparable-mass systems. E.g. in the equal-mass case (ly = i) the IGSF, ly- 
linear result (|2.17p . i.e. the contribution + 1.25 in Eq. dHI]), predicts Cl^f^ = 
6-3/2[1.3125] ~ 0.08930, which is significantly (8%) smaller than d'Cso'"^ = 
6-3/2[1.4294] ~ 0.09726, or the exact EOB value fjggf = 0.09670. A good 
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fit between EOB and equal-mass NR data requires that the equal-mass LSO 
frequency be roughly between 0.096 and 0.097. The approximate equation (|8.ip 
displays the needed complementarity between various approaches to the dynam- 
ics of binary systems: NR, GSF and EOB. 

4. We have also discussed (in Section 5) how the study of small-eccentricity 
orbits can allow one to confront the EOB formalism to GSF calculations. In 
particular, our Eq. (|5.2ip shows how the GSF computation, along the sequence 
of (quasi-)circular orbits, of the squared ratio between the radial (periastron to 
periastron) frequency and the azimuthal one gives us access to a function p{x) 
of X = (G(mi -I- 777.2) f^)^^'^ which, when interpreted within the EOB formalism, 
is an x-dependent combination of a(a;), a'{x)^ a"{x) and d{x). Here, a{u) = 
[d A{u]v)/dv\v=[) as above, and d{u) = [dD{u\v)/d 

T^jiy^O; where the EOB 

metric function I){r\ v) is defined in Eq. ()3.13p . More precisely, we recommend 
to study not only (as in [26]) the behaviour of p{x) near the LSO (where [ur/^Y 
vanishes), but also below the LSO (for ^ < a; < I), and around a; = 0. We 
have explicitly given the first two terms (2PN and 3PN) in the Taylor expansion 
(|5.32p of p{x) near x = 0. See Eq. I|5.39p (which also follows from results given 
in [39]). We emphasize the need of this comparison between 3PN results and 
GSF ones for confirming the validity of the results of [26| (similarly to the way 
3PN results were recently used [28] to confirm the validity of the gauge-invariant 
GSF result of [17j). We have also explicitly shown how the determination of 
the higher-order terms (p4 x* -t- x^ -f . . . , corresponding to 4PN + 5PN + . . .) 
in the Taylor expansion of p{x) can give us access to specific combinations of 
the higher PN coefficients 05, ae, . . . ; di, rfs, . . . entering the PN expansions of 
the two functions a{u) and rf(w), see Eq. (|5.33H5.37|l and (|5.45p . In particular if 
10 05 dominates over ^4, this could give us some interesting confirmation of the 
tentative determination (|4.46p of 05. 

5. We have also discussed several other ways to confront (conservative) 
GSF calculations to the EOB formalism. Some of them present challenges to 
the GSF line of work: such as the GSF determination of the total conserved 
energy and angular momentum (using either the second GSF approximation, or 
a detailed study of unbound orbits; see beginning of Section 6). An easier way of 
getting new, quantitative information from GSF studies is to study the special 
zero-binding zoom-whirl motion which starts, in the infinite past, with zero 
kinetic energy at infinite separation, and ends up, in the infinite future, whirling 
indefinitely around some limiting finite separation. We showed how the GSF 
study of this special motion could give us access to (at least) three dimensionless 
observables: j*, l^whiri and (a;^/il^)whiri- We related these observables to the 
values of various (i^-linearized) EOB functions at it = 1/4. We also ventured an 
approximate prediction for the value of a(|) entering the EOB prediction for 
j*, see Eqs. (jO^ and ([05)) . 

6. In addition, we have shown that the logarithmic terms In c that appear, 
at the 4PN level, in the post-Newtonian expansion of the near-zone metric of a 
gravitating system |44j give rise to a (Renormalization-Group-type) logarithmic 
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running of the expansion coefRcients of various functions: the "redshift" func- 
tion u*{y) of Ref. [27j (for which we computed the leading logarithmic term, 
Eq. (|7.14p ). and, arguably, the EOB function a{u) (and the other EOB functions, 
such as d{u)). We leave to future work an exploration of the effects of including 
such logarithmic terms, as in Eq. (|7.2p . in the EOB formalism. Though this 
will affect the details of our GSF/EOB comparison, we do not expect that this 
will introduce drastic changes because the EOB/NR comparison (and in par- 
icular the determination of the "good fit" region in the (05,05) plane) mainly 
depends on 05 and ag as effective parameters, describing the shape of the A{u) 
function in an interval between the LSO and the EOB light-ring (i.e. roughly 
between m ~ 1/6 and u ~ 1/3). Therefore, if we consider "running parameters", 
the EOB/NR comparison will be mainly sensitive to = (05 (In u)) and 

^effective _ (ag(ln u)), where the brackets denote an operation of averaging over 
an interval of u close to the LSO. As a rough approximation we expect that 
^effective ^ ^^^Yn ULSo) and agff'='="™ ~ a6(ln ulso)- Therefore, as the GSF re- 
sults depended on the behaviour of the a(u) function near the (unperturbed) 
LSO, the EOB/NR agff°'=*'™, af'"'*-^^'' should be also approximately relevant in 
the EOB/GSF comparison. 

Finally, let us remark that: 

a. We recommend that the GSF studies based on the use of the Lorenz 
gauge be systematically reformulated (or at least re-expressed) in terms of an 
"asymptotically flat" coordinate system. Indeed, the present use of a "non- 
asymptotically flat" coordinate system is not only confusing for general physi- 
cists, but can cause real errors (e.g. when considering inspiralling motions where 
the "renormalization" factor 1 -I- a, Eq. (|2.8p . connecting tLorcnz to iflat, would 
become adiabatically time-dependent). 

b. Our result (|4.46|) on the way the GSF result [26] breaks the degeneracy 
of the EOB-NR constraint on (05,05) implies, in particular, that the nonlinear 
dependence on v entailed by the Pade-resummed definition of the A(u\ v) func- 
tion within the EOB formalism [53] plays an important role. This can be seen, 
for instance, in considering an A function of the type 

^^-lincar^^. „■) = I _ 2u + IZ Op(u; O5, Og) , (8.2) 

where ap{u, 05, oe) is (uniquely) defined by requiring that the equal-mass value 
>lp''"°"(u; 1/4) be equal to the normal Fade- resumed ^pad|(^; ^5; ae) (whose 
structure was recalled in Section 4, before Eq. (|4.44p '). Then, one finds that, all 
along the banana- like "good fit" region in the (05,05) plane [M], the fractional 
i^-derivative (at = 0) cq, Eq. (|2.16p . of GMr^LSO predicted by the z^-linear A 
function (|8.2p stays in the vicinity of co [A''"''"''''''] ~ 0.82, practically indepen- 
dently of (05, og), when varying them along the center of the banana-like "good 
EOB/NR fit" region in the (05, og) plane. Note that the value 0.82 is only 65.6% 
of the GSF value (|2.17[) . By contrast, the value of cq corresponding to the exact 
(Fade) EOB function ^p2d|('"i ^5 '^5^ ^e) does vary along the central line of the 
"good fit" (05, Og) region, and does reach the GSF value (j2.17|) at (and only at) 
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the particular values (|4.46p . It would be interesting to confirm the need of such 
nonlinear behaviour in A{u; v) by exploring in more detail than has been done 
so far the EOB/NR comparison for several mass ratios. 
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